Unified description of floppy and rigid rotating Wigner molecules 

formed in quantum dots 
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Restoration of broken circular symmetry is used to explore the characteristics of the ground 
states and the excitation spectra of rotating Wigner molecules (RWM's) formed in two-dimensional 
parabolic A''-electron quantum dots. In high magnetic fields, the RWM's are fioppy rotors with the 
energies of the magic angular momentum (L) states obeying aL + b/L^^'^ . Under such fields the 
ground-state energies (referenced to the kinetic energy in the lowest Landau level) approach the 
electrostatic energy of A'^ point charges in the classical equilibrium molecular configuration. At zero 
field and strong interelectron repulsion, the RWM's behave like quasiclassical rigid rotors whose 
energies vary as L^. 

PACS numbers: 73.21.La; 71.45. Gm; 71.45.Lr 



The belief that the physics of two-dimensional (2D) 
semiconductor quantum dots (QD's) in high magnetic 
fields (B) is described 0, by composite-fermions, in- 
troduced 1^ originally for the bulk fractional quantum 
Hall effect (FQHE), has been recently challenged 4] by a 
proposal that the fundamental physical entity in strongly 
correlated QD's is a collectively rotating Wigner (or elec- 
tron) molecule (RWM or REM). Indeed, it has been 
lately found that electrons in 2D QD's can undergo in 
the strongly correlated regime a spontaneous phase tran- 
sition akin to the Wigner crystallization in the bulk, 
forming H, 0, S 13 specific polygonal geometric struc- 
tures that are called Wigner molecules (WM's). These 
geometric structures break the rotational symmetry (for 
symmetry restoration, see below) and reflect the presence 
of many-body crystalline correlations that are revealed, 
as illustrated previously, in the "humps" of the condi- 
tional probabilities ui Q and of the broken-symmetry 
electron densities IHUj- 

The majority of recent studies address the proper- 
ties of static WM's (as a function of the QD parame- 
ters) |l[iS[il[i3- 

Here, we study the properties of 
the rotating WM. We show through microscopic many- 
body investigations that the RWM in high i? is a floppy 
rotor, while at large Rw (and zero magnetic field) it 
transforms into a rigid rotor. [The Wigner parame- 
ter Rw = Q/hioo, where Q is the Coulomb interaction 
strength; Q = /{kIq), with Zq = ^?i/(m*a;o) being the 
spatial extent of the lowest single-electron wave function 
in the external parabolic confinement of frequency loq, 
and K is the dielectric constant.] The ability to capture 
the physics of the electrons in QD's in both the high 
magnetic field and the field-free regimes is an essential 
demonstration of the powerful unification offered by the 
RWM picture, a property not shared by other suggested 
approaches. 

The collective rotation of the WM is inherent and 
natural to the molecular picture. In particular, we 
show that it is manifested in characteristic energy-vs- 



angular-momentum relations for the excitation spectra 
and in specific limiting values for the ground-state ener- 
gies. These relations are important in the development 
of theories of electrons in QD's under the influence of 
a magnetic-field and in the field-free case, and they can 
be employed as diagnostic tools for assessing the validity 
and applicability of alternative theoretical descriptions. 

Description of the broad variation of the collective 
properties of electrons in 2D QD's requires a highly flexi- 
ble and accurate many-body method. In principle, exact 
diagonalization (EXD) could have been used; however, 
its computational limitations are a major obstacle. We 
demonstrate that the recently proposed 0, 0| two-step 
method of circular symmetry breaking at the unrestricted 
Hartree-Fock (UHF) level and of subsequent symmetry 
restoration via post-Hartree-Fock projection techniques 
is an accurate and computationally efficient approxima- 
tion which provides a unified microscopic description of 
the emergent picture of RWM's in 2D QD's. 

The restoration of broken circular symmetry is neces- 
sary for a proper description of RWM's. Indeed, while the 
broken symmetry UHF solutions describe static WM's 
[3, llfl llll I it is the symmetry restoration step which 
describes the rotation of the WM's and underlies the dif- 
ferentiation between rigid and floppy (at high B) rotors. 

The two-step method. In general, the localized broken 
symmetry orbitals are determined numerically via a self- 
consistent solution of the UHF equations . An efficient 
alternative, however, is to approximate these orbitals by 
appropriate analytical expressions 0, 0, 0| . Since 
we focus here on the second step (restoration of the circu- 
lar symmetry) and the case of high B, it will be sufficient 
to approximate the UHF orbitals (first step of our proce- 
dure) by (parameter free) displaced Gaussian functions 
lll^ : namely, for an electron localized at Zj, we use the 
orbital 
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FIG. 1: Total energies, E — 2hi}, for the triplet state of 
N=2 electrons as a function of the magnetic field B. Solid 
line: PRJ (rotating WM); dashed line: exact; dotted line: 
UHF (static WM). The integers below the troughs of the PRJ 
and exact lines indicate the corresponding magic angular mo- 
menta. The parameters are: m* = 0.067me, Tiuq — 3 meV, 
and K = 12.9. Here, as well as in the other figures, the Zeeman 
contribution is not included, but can be easily added. 
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FIG. 2: PRJ (solid lines) and UHF (dashed hues) total en- 
ergies, E — NhQ,, for the fully spin-polarized ground states of 
A'^ = 3 and N — 5 electrons, as a function of 1/B. The range 
covered is 4 T < _B < 400 T, and the zero of energy corre- 
sponds to Ell = 14.415 meV for iV = 3 and to E^l = 42.873 
meV for N = 5. The integers below the troughs of the PRJ 
curves indicate the corresponding magic angular momenta. 
The parameters are as in Fig. 1. 



with z = X + iy, Zj — Xj + lYj, and A — ^/h/rn*Q; 
O = y/ujQ + oJc/i, where uJc = eB/ {m*c) is the cyclotron 
frequency. The phase guarantees gauge invariance in the 
presence of a perpendicular magnetic field and is given in 
the symmetric gauge by (p(z, Zj; B) = [xYj — yXj)/2l'^, 
with Ib = \/hc/eB being the magnetic length. We only 
consider the case of fully polarized electrons, which is 
appropriate at high B. 

We take the Zj's to coincide with the equilibrium 
positions (forming nested regular polygons ^14]) of N 
classical point charges inside an external parabolic con- 
finement of frequency cuq, and proceed to construct the 
UHF determinant *UHF|-^j orbitals u{zi, ZjYs, 

i,j = 1,...,N. Correlated many-body states with good 
total amgular momenta L can be extracted from the 
UHF determinant using projection operators; the pro- 
jected energies are given by |6| 

EpKj{L)^J^ M7)e*T^d7y^ n(7)e'^^d7, (2) 

with hi-f) = («'UHF(0)|iJ|5-UHF(^)^ ^ 

(^fUHF(o)l^UHF(^)^^ where ^uhf^^) original UHF 

determinant rotated by an azimuthal angle 7 and H is 
the many body Hamiltonian (including the vector poten- 
tial, external confinement, and Coulomb two-body repul- 
sion). We note that the UHF energies are simply given 
by £;uHF = h{0)/n{0). 

Compared to an EXD calculation, the CPU time for 
calculating the projected energies [see Eq. (0)] increases 
much slower as a function of L and N. In addition, the 
computational efficiency of Eq. Q is greatly enhanced 
by the knowledge ^3 of the analytical forms of the (in 
general complex) matrix elements of the single-particle 



and two-body components of H between the displaced 
Gaussians. 

To test the accuracy of our method, we display in Fig. 
1 the total energies E — 2h^ as a function of the magnetic 
field for the triplet state of A'' = 2 electrons and for all 
three levels of calculations, i.e., the UHF, the subsequent 
projection [PRJ, see Eq. and the exact result f^. We 
observe that the UHF energies fare poorly compared to 
the PRJ ones; note that for _B > 7 T, the PRJ and exact 
results are practically the same. We further note that the 
UHF energies do not reproduce the characteristic oscilla- 
tions of the exact ones; these oscillations originate from 
the fact that only states with magic angular momenta 
can become ground states (in the case of the triplet state 
for N — 2, the magic angular momenta are L = 2m + 1, 
m = 0,l,2,...). 

Ground states at high magnetic field. Fig. 2 displays 
the ground state energies E — Nhi^ for = 3 and N = 5 
electrons as a function of 1/B. Since our method allows 
us to reach magnetic field values as high as i? = 400 
T, we conclude unequivocally that, as -B — s- 00, both the 
PRJ and UHF energies approach the limiting value corre- 
sponding to the classical energy of N point-like electrons 
at equilibrium [in the (0,N) configuration for = 2 — 5] 
inside an external parabolic confinement of frequency loq, 
i.e., 

KlW = (3/8)(2i?w)'/'^5^/'?iwo, (3) 

with Sn = {sin[(j - l)7r/7V])"^ d- Except for 

the B ^ 00 limit, the UHF energies are higher in value 
than the projected ones. 

The projected ground states have magic angular mo- 
menta that vary (from one trough to the other) by steps 
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FIG. 3: Projected total energies, iSpRj — 5hQ,, for the yrast 
rotational band (states with the lowest energy at magic an- 
gular momenta) of A'^ = 5 fully-polarized electrons at -B = 60 
T, as a function of the angular momenta L. The ground state 
has angular momentum L = 115. The choice of parameters is: 
m* = 0.067me, hu>o = 4 meV, and k — 12.9. Inset: The two 
contributions i5pj^j(I/) — 5hQ, and Vprj(Z/) associated with 
the confinement (lower curve) and the Coulomb interaction 
(upper curve), respectively. 



of N units (A^ — 3,4,5) in accordance with the trend 
found for a much smaller range of angular momenta 
in earlier EXD studies d, 0|. With our projection 
method, we were able to reach remarkably high values 
of L = 369,690, and 1100 for iV = 3,4, and 5, respec- 
tively. The fact that E'prj — Nhfl and inexact — Nhfl 
(see Ref. 0) tend to Ell{N) for B — !■ CX3 is a significant 
finding that has not appeared in the FQHE literature 
{L JJ.\ ; instead it has been implied that E — NhQ ^ as 
B oo. In any case, it is doubtful that the composite- 
fermion approach can conform to this rigorous limit 'l6"| ; 
indeed the precise value of E^l{N) depends crucially on 
the long-range character of the Coulomb force acting be- 
tween localized classical particles, and thus it is a result 
innate to the WM picture. 

Furthermore, this limit points out to interesting 
experimental ramifications. Namely, at high B, 
the addition energies (expressed by the finite sec- 
ond difference A2£'(A^)), while not conforming to the 
capacitance— Coulomb blockade model will tend to 

r,^l{N+l) + E^l{N-l)~ 
3.446 meV for = 4 and the parameters 
listed in the caption of Fig. 2); this constant value maybe 
detected experimentally. 

The yrast rotational band at high B. The yrast band 
consists of the lowest energy states (for a given B) having 
magic angular momenta (the ground state is a member 
of this band and the remaining states are excited ones). 
The yrast band for = 5 calculated at i? = 60 T is 
displayed in Fig. 3. 

The energy of the yrast states is composed Ij of two 



a constant value /!^2c\^{^) 
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TABLE I: Projected interaction energies Vprj(I/) associated 
with yrast states at two values of B = 60 T and _B — » oo 
for A^ = 5 [the (0,5) configuration] and A^ = 6 [the (1,5) 
configuration], respectively. / = Vprj{L)L^^^ /Cv, where 
Cv = 15.388 for Af = 5 and Cv = 18.787 for A^ = 6. The 
fractional fillings are calculated through f = N{N — 1)/{2L). 
Energies in units of c^/kX and c^/kIb for _B = 60 T and 
B — > oo, respectively. At B ^ oo, V{L) is independent of 
hiJo- At large finite B, V{L) depends on hu>o through A. 
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contributions; namely, E{L) = E'^{L) + V{L), where 
E'^'{L) is the confinement energy associated with the ex- 
ternal potential (including the kinetic energy), and V{L) 
is the interaction energy due to the Coulomb force. For 
high B, it is well known P, ^3 that the confinement 
energy varies linearly with L; in particular, E'^{L) — > 

n{n-ujj2)L + Nnn. 

We have calculated the projected values [see Eq. 
^] for these two components, namely, E^^j{L) and 
VpRj(L), for = 2 — 6 and a variety of B values 
and QD parameters. Our ability to calculate efficiently 
high L values allowed us to confirm numerically the lin- 
ear L-dependence of the confinement energy E'^{L), and 
also the precise proportionality coefficient 7i(r2 — ujc/i). 
Moreover, rather unexpectedly, we also find that the in- 
teraction energy behaves asymptotically as L~^/^, i.e., 
V{L) Cv/L^^^: where Cy is a constant given below. 

TABLE I lists Vprj(L) for TV = 5 and B = 60 T, as 
well as for = 6 [the (1,5) configuration] and B 
oo; the latter limit can be easily taken by setting 0] 
A — IbV2 in Eq Q, which restricts the electrons to the 
lowest Landau level. From TABLE I, one can see that the 
quantity / = VpRj(L)i^/^/C\/ approaches unity as the 
angular momentum intreases, where Cy — K^^^{Sk/4^ + 
Sk^n-i) e^/inX), with K = N or N - 1 iov WM's in the 
(0, iV) or (I,iV — 1) configurations, respectively. 

The Cv/L^/"^ asymptotic dependence of the interac- 
tion energy can be readily associated with a model Q 
of a classical floppy molecule rotating inside a parabolic 
confinement characterized by a frequency Q, [l^ . Indeed, 
the energy of such a molecule is given by 
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TABLE II: Projected total energies _EpRj(L) at B = and 
Rw ~ 200 associated with yrast states for N — 5 electrons 
[(0,5) configuration]. / = Cr/Cj^ (see text). Energies in 
units of fiLoo. 















L 


-EpR.7 


/ 


L 


-EpRJ 


/ 





323.3070 




25 


324.7657 


0.988 


5 


323.3656 


0.992 


30 


325.4033 


0.986 


10 


323.5414 


0.992 


35 


326.1537 


0.983 


15 


323.8338 


0.991 


40 


327.0153 


0.981 


20 


324.2422 
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where J {a) = Km*a^ is the moment of inertia of K 
point-like electrons located at the vertices of a regular 
polygon of radius a. The second term is the potential en- 
ergy due to the confinement, and the last term is the clas- 
sical Coulomb-repulsion electrostatic energy. At given B, 
the radius of this floppy molecule varies with L, and for 
large angular momentum (and/or high B) it is given by 
a « \y/L/K [a ^2ncL/{eBK) for B -> oo]. Substi- 
tution into Eq. (0J yields the aforementioned linear and 
contributions, i.e., 

Elf{K) « n{Vl - Luj2)L + (5) 

A semiclassical approximation Lg' of the ground-state 
angular momentum minimizes Eq. (O, yielding oc B. 
Consequently, the ground-state radius ags as B ^ 

oo, where rg = IqR}^{Sk/^ + 5k,n-i)^^'^ is the radius 
of the static classical Wigner molecule. Note that, in the 
ground state as B increases, the WM rotates faster in 
order to maintain the requisite value of ags = tq. 

We note that this inverse-square-root-of-L law at con- 
stant B for the Coulomb interaction energy has also been 
overlooked in the FQHE literature 1, 17]. This law pro- 
vides a nontrivial test in that it lends further support to 
the RWM picture; especially since it also applies in the 
B ^ oo limit (see right part of TABLE I), which restricts 
the electrons to the lowest Landau level (FQHE regime) . 
We note that high L's correspond to lower fractional fill- 
ings, experimentally achieveable in high- mobility sam- 
ples. 

The rigid rotor at zero magnetic field. At _B = 0, we 
found it advantageous to allow the width A and the po- 
sitions Zj 's of the displaced Gaussians to vary in order 
to minimize (for each L) the projected energy [Eq. 0]. 
For large Rw^s and for TV = 2 — 5 electrons, we find 
that the energies of the states in the yrast band can be 
approximated by 

SpRj(L)«SpRj(0) + CflL2^ (6) 

where the rigid-rotor coefRcient Cn is essentially a con- 
stant whose value is very close (see the / values) to 
the classical one, corresponding to point charges in their 
(0, TV) equilibrium configuration inside a parabolic con- 
finement of frequency luq] i.e., Cr ~ Cj^ = h? /{2J{ro)). 



In Table ^ we list calculated i5pRj(L) values for 
N = b when Rw = 200. The fact that / = CR/Cj^ « 1 
for L < 45 illustrates that the RWM behaves as a quasi- 
classical rigid rotor. For smaller values of Rw, the rigid- 
ity of the RWM is progressively reduced . We note that 
-E-PRJ (0) — NhujQ is very close to the classical electrostatic 
value given by Eq. (jSJ. 

Conclusions. Using the method of restoration of bro- 
ken circular symmetry, we have shown that the RWM's 
in QD's exhibit characteristic properties as a result of 
their collective rotation: (I) in high B, they behave like 
floppy rotors and the energies of the magic states in the 
lowest rotational (yrast) band, associated with magic 
angular momenta, obey aL + b/L^/"^. In addition, the 
ground-state energies (referenced to the kinetic energy 
in the lowest Landau level) approach the limit of the 
electrostatic energy of N point charges in their classical 
equilibrium molecular configuration; (II) at i? = and 
strong interelectron repulsion, the RWM's behave like 
quasiclassical rigid rotors whose yrast band exhibits an 
L^-dependence. The rotating- Wigner-molecule picture 
unifies the description of the various physical regimes of 
strongly correlated electrons in 2D QD's. This unifica- 
tion is achieved through the computationally powerful 
method of restoration of broken symmetries via projec- 
tion techniques. 

This research is supported by the U.S. D.O.E. (Grant 
No. FG05-86ER45234). 
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